%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Calcul de l'interpolation des moments la poutre
%   entrée:

%   sortie:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [momentInt] = Moments_interpoles (donnees)

momentInt = [];
momentIntTemp = zeros(1,donnees.nInt+1);

figure;
hold on;
for i=1:(length(donnees.maillage)-1)
    
    elemTemp = donnees.maillage(i:i+1);
    pas = (donnees.maillage(i+1)-donnees.maillage(i))/donnees.nInt;
    UhElem = donnees.res.Uh(2*i-1:2*i+2) ;
    
%     for j=0:(donnees.nInt)
%         
%         momentIntTemp(j+1) = donnees.poutre.E*donnees.poutre.I*(...
%             donnees.res.Uh(2*i-1)*phiISec(j*pas,elemTemp) +...
%             donnees.res.Uh(2*i)*psiISec(j*pas,elemTemp) +...
%             donnees.res.Uh(2*i+1)*phiJSec(j*pas,elemTemp) +...
%             donnees.res.Uh(2*i+2)*psiJSec(j*pas,elemTemp));
%     end

    for j=0:(donnees.nInt)
        
        momentIntTemp(j+1) = Moment_element(donnees,j*pas,elemTemp,UhElem) ;
        
    end
    momentInt = [momentInt;momentIntTemp];
plot( donnees.maillageInt((i-1)*donnees.nInt+1:i*donnees.nInt+1) , momentIntTemp )   
end

end

